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Abstract 

We construct and investigate a family of two-band unitary systems living 
on a cylinder geometry and presenting localized edge states. Using the trans¬ 
fer matrix formalism, we solve and investigate in details such states in the 
thermodynamic limit. Analitycity considerations then suggest the construc¬ 
tion of a family of Riemman surfaces associated to the band structure of the 
system. In this picture, the corresponding edge states naturally wind around 
non contractile loops, defining by the way a topological invariant associated 
to each gap of the system. 


1 Introduction 

First discovered in the context of the quantum Hall effect boundary states 

turn out to be the hallmark of topological properties that can emerge in any di¬ 
mensions and for various symmetry classes HE]. It was also realized that these 
topological properties can be achieved in miscellaneous physical systems beyond 
solids PE] leading for instance to the discovery of chiral edge states in cold-atom 
(HI, electromagnetic P ITTH ITT] and acoustic Ha [131 n E5] lattices. 

These topological properties characterize the bulk bands and are encoded by 
a topological invariant whose integer value cannot change unless bands touch, or 
equivalently, unless the gap closes. For instance, in two dimensions and in absence 
of particular symmetry, this topological invariant is the hrst Chern number uniiiz]. 
As the edge states live in the gaps spectrum, it follows that they cannot be removed 
or added unless a topological transition of the bulk bands happens when the gap 
closes. This gives a topological robustness to the edge states. 

Recently, similar topological properties have been found in unitary systems, 
namely physical systems whose behavior is described by a unitary operator rather 
than a Hermitian operator. It follows that their spectrum is periodic, unlike an 
energy spectrum which is bounded. Among them are the Floquet systems that 
are periodic in time [T8l 119] . These dynamical systems, such like periodically driven 
solids |2ni[2T], shaken cold-atom gases |22j, photonic lattices j23l[2l| or discrete-time 
quantum walks (23] are fully described by their (unitary) time evolution operator. 
Importantly, it was shown that beyond the analogy with equilibrium (Hermitian) 
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systems, they can exhibit edge states whereas the usual bulk topological invariants 
vanish. This is understood as an exotic topological property characterized by a 
novel invariant assigned to a gap (rather than a band) and that accounts for the 
time evolution over a period 123 EH- 

However, there are unitary systems which are not Floquet systems but still ex¬ 
hibit edge states. This is the case of photonic and microwave networks which can 
instead be described by (unitary) scattering matrices [28l |29l 130]. How to describe 
the topological origin of these edge states? Can it still be related to a bulk property? 

To answer these questions we construct, in section]^ a generic model of a two-gap 
unitary system on a cylinder geometry. This model achieves the specihc situation for 
which the topological invariants of the bands (namely the Chern number in our case) 
vanish. Then, in section we apply the transfer matrix method to investigate the 
appearance of edge states and discuss their simultaneous existence in the two gaps 
(thus guarantying a vanishing Chern number of each band). In its seminal paper, 
Hatsugai showed that the transfer matrix approach provides a deep understanding 
of the topological nature of the edge states in the quantum Hall phase |31j . In 
particular, the transfer matrix allows one to focus directly on the gaps where the 
edge states live rather than the bands only. We follow the same strategy and adapt 
this method to the unitary case in section Our analysis reveals the underlying 
topologically non-trivial structure of the edge states and justihes the dehnition of a 
topological invariant assigned to a gap instead of a band. Finally, several examples 
of physical systems ruled by the present model are discussed in section 


2 Two-gap unitary models with topological edge states 

2.1 Heuristic construction on a finite size lattice 


We construct heuristically a two-gap unitary model for a strip geometry that exhibits 
topologically protected edge states. For simplicity, we shall treat the cases of 0 or 1 
edge state, but the generalization to several edge states is straightforward. To this 
end, we consider a system with two degrees of freedom - that we refer to as A and B 
- in the cylinder geometry such that the (dimensionless) quasi-momentum k G U{1) 
is a well-dehned continuous parameter in the periodic direction, whereas the lattice 
remains hnite in the other one, as sketched in hgure (a). A state |d'(A:)) of the 
system is then given by a 2Wcomponent vector (y4i(/c), Bi{k ),..., A^l^k), B]s[{k))'^, 
which, for the scope of this study, is an eigenvector of a unitary matrix U{k) G 
U{2N), that is 

U{k)\^{k)) =e-^^^^^\^{k)) . (1) 

We would like the phases e{k) G Si of the eigenvalues of U{2N) to display two-gaps 
and two edge states lying in these gaps. Formally, the simplest way to obtain such 
a phase spectrum is to impose a quasi-diagonal form for U{k) such as 


/eifc 

Uoik) = lay® In-1 


( 2 ) 
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where the bulk part iay yields two flat bands e = ±7r/2 (in the range e g] — tt, tt]), 
whereas e**^ and e“*^ guaranty the existence of propagating modes at the boundaries 
Ai and respectively, (see figure (b)). The phase spectrum e{k) of Uo{k) can 



Figure 1: (a) Sketch of the cylindrical geometry for the strip of width N. The two arrows 
represent the chiral edge states at the boundary n = 1 (purple) and N (red). Phase 
spectrum e(k) obtained by a direct diagonalization of U{k) with = 10 for 0 = 0 (b) and 
0 = vr/S (c). 


then be obtained as a particular case of a more general unitary matrix 


t' d 

'a Hi 


U{k)=\ U2 ®In-i 

where G U (2) and 1 is a scalar. One gets explicitly 


/ 


U{k) = 


d P'l 




J 


V 


, A = 


Pld -rid t[p'^ p[p'^ 


P1P2 rip 2 t[T2 P\T2 


Pi n y 

where, for concreteness, we choose 

t[ = tI = cos 9 , di = ~Pi = ^2 = T 2 = sin 9 , P 2 = ~P 2 = cos 9 . 


(3) 


(4) 


(5) 


In Eq. (Q, A is a 2 X 4 matrix coupling components A and B at sites n and n + 1 for 
1 < 77 , < (hence describing the bulk properties), whereas first and last lines are 
constraints for A and B at sites 1 and N (hence yielding the boundary conditions). 

This model allows us to study the fate of the edge states when varying 9. First, 
note that lJ{k) corresponds to f/o(A:) for 0 = 0 , up to a unitary transformation 
that does not change the spectrum. Then, for 0 7 ^ 0, the bulk bands acquire a 
dispersion (figure (c)) and eventually touch at (A; = 0, e = 0) and (A; = vr, e = vr) 
for 0 = 7 r/ 4 , thus implying a transition toward a gapped phase with no edge state 
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( 7 r /4 < 9 < 7 r/ 2 ). A similar matrix U{k) was derived to investigate the propagation 
of electromagnetic modes in arrays of optical resonators 12 ^ . and can also me 
adapted to describe discrete-time quantum walks for a spin- 1/2 particle or either a 
time-dependent tight-binding model, as we discuss in section 

In the following, we investigate the topological properties of this model with the 
transfer matrix formalism. 


2.2 Transfer matrix formalism 


The transfer matrix formalism is a standard method to tackle a large variety of 
problems. The starting point consists in relating the wave function amplitudes on 
adjacent sites of a lattice by a matrix T. In our case, this translates as 


\ -Dn+l / \ -Dn / 


( 6 ) 


From the general form of the matrix U{k) in Eq. Q, one can infer the relation 


AP2 

t'iT 2 - e~ 


PlP2 


An+i 

Bn+l 


-P 1 T 2 -TiT^ e“ 
—Plp2 —'T 1 P 2 


B„ 


(7) 


Note that the relation Q, valid for 1 < n < describes the bulk; it does not take 


into account the boundary conditions (n = 1, N) that will be given in section 3.2.1 


We can then deduce the expression of T which, for the specihc parameters (|^, reads 

tan^e''^ —pk^-fe) 


Te{Ke) = 


—ekc+fc) _|_ tan6' 


cotan 6 ' e^*^ -|- 


e*^" ""^-l-tan^ 

—a — 2 cos k 

sm 6 cos 6 


( 8 ) 


By construction, the matrix T* describes how the amplitude of a state evolves 
from the edge n = 1 to a site n = I of the bulk. This information is then encoded into 
the eigenvalues A±(e) of the transfer matrix T. In particular, either |A±(fc,e)| = 1, 
which corresponds to an oscillatory (or delocalized) bulk state at {k,e)] or one of 
the two eigenvalues satishes |A(/c,e)| < 1 - note that A+ = Al^ since detT = 1 - 
and no delocalized state can appear i.e. there is a band gap. The eigenvalues of T 
being solutions of the characteristic polynomial A^ — tr TA -|- det T = 0 , they read 


f>l: \i(k,e) = -f±^P-l, 

A±(*:,t) = -/±i\/l-/V 

where / = fe{k,e) = —\tiT{k,e) that is 

f = cos k- — -- . 

sm 9 cos 9 


|A±(A:,e)| 7 ^ 1 (9) 

|A±(fc, 6 )| = l (10) 

( 11 ) 


In other words, for the model ([^, the eigenvalues of the transfer matrix are fully 
determined by its trace. Finally, the two eigenvectors v__ and v+, associated re¬ 
spectively to the eigenvalues A_ and A+ of the matrix T = {T/j}, can be written 
as 

^ (a± -Vii) • 

In the following, we show how the projected bulk bands, the edge states and 
their topological properties can be inferred from these eigenvalues and eigenvectors. 
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3 Band structure from the transfer matrix 


The transfer matrix formalism allows one to solve explicitly the initial problem Q, 
namely to reconstrnct the band strnctnre and the corresponding states. We focns 
on the thermodynamic limit —)■ oo of the infinite cylinder: in that case the bnlk 
bands become a continnnm, whereas the edge states persist. 


3.1 Bands and gaps 

The bnlk bands e{k) correspond to regions of the diagram [k, e) for which the eigen- 
valnes of the transfer matrix satisfy |A±(fc,e)| = 1. For the model we consider, this 
region is determined by the solntions of f‘^{k,e) ^ 1 according to Eq. (10). This 
can be fonnd explicitly by using the expression of f{k, e) in Eq. (0. The result is 
shown in figure and compared to a direct diagonalization of U. 




k 


Figure 2; Bulk bands of the strip obtained (left) from a direct diagonalization of U{k) 
with = 10 and (right) from the transfer matrix. The gap/band edges dispersions are 
represented in blue (red) for the gap g- {g+). The shaded area corresponds to the bulk 
band continuum at N —)> oo. 


Consistently, the gaps correspond to regions of the diagram (fc, e) for which the 
eigenvalues of the transfer matrix satisfy |A±(fc,e)| ^ 1, that is when p{k,e) > 
1. This inequality has two solutions f{k,e) < —1 and f{k,e) > 1 yielding two 
different domains in the diagram {k,e). Since these two domains are separated by 
a band region (/^ ^ 1), they thus correspond to the two gaped regions that we will 
refer to as gap when / < — 1 and when / > 1. The gaps g^ and g+ are 
respectively delimited by/ = —lor/ = +l. We therefore end up with a criteria 
that distinguishes the bands and the two gaps directly from the trace of the transfer 
matrix 


/ < —1 -v^ gap g- , ^ 1 -v^ bulk bands , / > 1 gap g+ (13) 


The explicit dispersion relation of the two gap edges are then inferred from Eq. (11) 


G-^t{k) = arccos (sin6'cos6'(cosfc + 1)) 
G+,b(fc) = arccos (sin 6 cos 6'(cos k — 1)) 


G-^b{k) = — arccos (sin 9 cos ^(cos k + 1)) 
Gj^ t{k) = — arccos (sin 6^ cos 6*(cos k — 1)). 

(14) 
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and plotted in Fig. 


3.2 Edge states 

3.2.1 Boundary conditions 

So far, the bulk part of U have been encoded within the transfer matrix. To charac¬ 
terize the existence of edge states, one also need to specify the boundary conditions. 
Whereas in Hermitian systems, a standard boundary condition on a lattice consists 
in imposing the vanishing of the wavefunction, such a procedure does not apply here 
since it is not compatible with unitarity. 

As for the bulk transfer matrix, the boundary conditions are instead inferred 
from the original eigenvalues problem on U. From Eq. (|^, we relate for each edge 
n = 1,N the amplitudes An and Bn with the phase e and the coefficients of Ui as 


{t[ - e Ai = -p'lBi , piAn = (e - ri) Bjy . (15) 


A geometric way to reformulate these boundary conditions is to introduce two bound¬ 
ary vectors bi and that are parallel to ai = {Ai,Bi) and Un = that 

is 

det(ai,bi) = 0 , det(aN,bN) = 0 (16) 

where we have defined 


bi = 





(17) 


These two boundary vectors are related by the transfer matrix as T^“^bi oc bN- 
This equation constrains the allowed values of couple (fc, e) for finite size A^, thus 
yielding the corresponding solution of the initial problemj^ ([^ . Next, we impose that 
in the thermodynamic limit (iV —)■ cx)), bN does not depend on N anymore. Thus, 
in that limit, we can write T^bi = so that bi becomes an eigenvector of 

the transfer matrix, given by (12). Then, by decomposing bi as bi = av_|_ /3v_, 

we end up with 


(A - 1) v+ + (A-^ - l) v_ = 0 (18) 

(where we have set A+ = A for commodity). Let us assume in addition that |A| < 1. 
Then |A|^“^ vanishes and tends to infinity. It follows that /3 = 0 and thus 

bi is proportional to v_|_. We can show in the same way that bN is proportional to 
v_ when N ^ oo. This result will allow us to fully characterize the edge states in 
the next section. 

Besides, ai and Un beeing also eigenstate of T in that case, the corresponding 
state |T) of the initial problem is exponentially decreasing or increasing with typical 
length In A, namely it is localized at one edge of the system. 

^In particular this equation also gives the bulk bands dispersion relation at finite size. 
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3.2.2 Existence and dispersion relation 

Concretely, in order to know which one of the A±’s satisfies |A| < 1, one has to fix /, 
which amounts to fix the gap. Let us for instance focus on the gap (i.e. / < — 1). 
Then, it is clear that |A+| > 1 and thus |A_| < 1. Therefore, in the gap g^, an 
edge state at n = 1 is given by the proportionality relation between the boundary 
vector bi and the eigenvector v_, whereas the edge state at n = is given by the 
proportionality relation between the boundary vector Ln and the eigenvector v+. 
The opposite situation occurs in the gap g^, for which / > 1, and thus |A_| > 1 and 
|A+| < 1, meaning that an edge state localized at the edge n = 1 (n = A^) is now 
associated with the eigenvalue A+ (A_) and the eigenvector v+ (v_). The existence 
of an each edge state and its dispersion relation is therefore given independently for 
each gap and for each boundary n = 1, A^ by the zeros of the function 


In,±{k,e) = det(b„, v±) 


(19) 


as summarized in table [H 



gap g- 

gap g+ 

At edge n = 1 

Ii-{k,e) = 0 

h,+ {kA) = 0 

At edge n = N 

lN,+ ik,e) = 0 

lN-{k,e) = 0 


Table 1; Equations for the existence of edge states in the space parameters {k, e). 

Let us treat the case of the edge state n = 1. Each of the two equations 
h,±{k,^) = 0 can be split into two equations corresponding to the real part and 
the imaginary part of e). These two equations can be treated simultaneously, 

and after some algebra, one gets 


Q{k,e) = 


sin e = — cos 6 sin k 
cos e cos k 


cos 9 sin 9 


= ±y7rrT, 


( 20 ) 

( 21 ) 


Whereas the explicit dispersion relation of the edge state localized at edge n = 
1 is straightforwardly obtained by inverting Eq. (20), its domain of validity in 


k is constrained by the equation (21) which depends on the gap g± through the 


contribution ± of the square root. Concretely, the edge state only exists in the 
region (fc,e) satisfying by Q{k,e) > 0 for the gap g+ and in the region Q{k,e) < 0 
for the gap g^. A similar procedure can be performed for the other edge so that 
the spectrum of the strip in the thermodynamic limit is recovered. As displayed in 


figure 

figure 


3 this procedure provides the dispersion relation of the edge states (see also 
6 (a) for the reconstruction of the full spectrum). Moreover, it shows that 
edge states cannot exist if the sign of Q does not change in the gap. It reveals the 
importance of the role of the branch ± of the square root in Eq. (21). This is a 


crucial point to understand the topological robustness of the edge states that will 
be developed in section 
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Figure 3: Dispersion relations for (a) 9 = Stt/S, (b) 9 = 7r/4 and (c) 9 = vr/S. Bine area: 
Q > 0 in the gap, pink area: Q < 0 in the gap and dashed line: Q = 0. The edge state 
(represented here at bonndary n = 1) exists and tonch a band at k = ± arccos (± tan^ 0) 
when Q changes sign. 


3.2.3 Gaps correspondence 


Before discussing the topological nature of the edge states, we emphasize that they 
always appear simultaneously in the two gaps for any generic unitary matrix U of 
the form of Eq. ([^ (provided the spectrum of U is gapped). This is a particularly 
interesting property of unitary systems since the Chern number of the Bloch bands 
of the periodic system in both x and y directions is guaranteed to vanish in this case 
|26] . However, the conditions to achieve such a remarkable topological property are 
not clear and remain an open question. Here, we point out that the transfer matrix 
formalism is particularly useful to reveal a correspondence between the gaps that 
constrains the existence of an edge state simultaneously in each gap. 

To make explicit such a property, we define a map v acting in parameters space 
V : {k, e) —)■ (nfc, ve) that relates an edge state at the boundary n = 1 that exists in 
the gap g+ to an edge state localized at the same boundary but existing in the gap 
g^. We aim at determining under which conditions e) and In,-{vk, ve) vanish 

simultaneously. By swapping the gaps, the function / changes sign as f{vk,ve) = 
—f{k, e) so that the eigenvalues become 

X±{k,e) —>■ X±{vk,ve) = —A=p(fc,e) . (22) 


It is clear that a map v acting on the eigenvalues of the transfer matrix as in Eq. 
(22) actually changes the sign of the trace of the transfer matrix. This allows us to 
express such a constraint on the transfer matrix itself as 


o'z T{k, e) Gz = —T T{vk, ve)T 


(23) 


with Gz the standard Pauli matrix and where T is allowed to be either the complex 
conjugation operator k or the identity. It is then easy to check that, for the model 
([^, the map n : (fc, e) —?• (vr — fc, tt — e) (defined modulo 27r) satisfies Eq. (23) with 
T = K, so that Eq. (22) is satisfied as well. It follows that an eigenvector v±(/c,e) 
of the transfer matrix T{k,e) (see Eq. (12)) is transformed as 

v±(nfc,ne) = GzK\zp{k,e) (24) 
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since Ti 2 {vk,ve) = 7^2(^!^) A±(fc,e) = A^(A:,e) for the model (|^. 

A relation of proportionality between /i ±(/c,e) and Ii^z^{vk,ve) can finally be 
inferred provided that the boundary vector bi(fc,e) - being itself an eigenvector of 
the transfer matrix in the thermodynamic limit - is also transformed as 


h-i^{vk,ve) = (T^Kbi(A;,e) (25) 

(as it can be checked explicitly from Eqs. ([^ and (12)). Indeed, according to Eqs. 
and ([25| one gets 


h,±{k, e) = K det (a^bi, cx^Vzp) {vk, ve) (26) 

= —K det (bi, V::p) (t;fc, r;e) (27) 

=—K Ii^^p(vk, ve) . (28) 

and thus /i,±(/c,e) and Ii^zp(vk, ve) vanish simultaneously. This shows that if an 
edge state localized at edge n = 1 exists in the gap g+ (defined by the zeros of /i,+), 
then so does an edge state localized at the same edge in the gap (dehned by the 
zeros of A _). 


4 Topological property of the edge states from the 
transfer matrix 

4.1 Riemann surface 

4.1.1 Motivations 

We would like to specify the topological nature of the edge states regardless the 
topological properties of the bands which are anyway always characterized by a 
vanishing Chern number as discussed above. A natural approach then consists in 
assigning a winding number to the edge states themselves [3^ |3Sl El]- When 0 = 0, 
then U{k) reduces to a block-diagonal matrix where the contributions of the bulk 
and of each edge can be distinguished. In that case, we notice that each edge term 
h±{k) = e^'^ is a 7/(1)—valued scalar so that its winding is well dehned [26] and we 
get ^ dk h}^dkh± = ±1. How to construct such a quantity for arbitrary 9 when 
U{k) is not block-diagonal anymore? How to assign a specihc gap to this winding? 
It is clear that the dispersion relations of the edge states obtained before will be a 
key ingredient. However they are not dehned for every k, not even periodic when 
considering one edge, and multivalued when considering both edges in one relation. 
Hence the winding number of such quantity is ill dehned. 

A possibility to overcome this problem is to consider instead the functions 
In,±{k,e), whose zeros fully dehne the edge states dispersion relations. These func¬ 
tions depend explicitly on the eigenvalues A± of the transfer matrix that contain the 
square root a//^ — 1. Therefore, their hrst derivative diverge around / = ±1, that 
is when [k, e) touch the bands edges. This divergence of the hrst derivative appears 
in the wave function too, which - remembering that the boundary vectors are also 
eigenvectors of the transfer matrix - is proportional to (An, Bn)'^ oc A(^“^bi_Ar. Thus, 
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to properly define a topological invariant, one first needs to find a way to obtain 
smooth solutions for the wave-functions when changing the branch of the square 
root ± a //2 - 1 = Fa //2 - 1. 

The natural tool to deal with this issue is the theory of Riemann surfaces, as 
used by Hatsugai in 1993 to correctly define the winding number of the edge states 
in the (Hermitian) Harper model j3l]. To construct the Riemann surfaces associated 
to our problem, we adapt the approach by Hatsugai to the unitary case where the 
eigenenergies are replaced by exp(—ie) which are t/(l)-valued. Substituting / by its 
expression, this approach suggests to consider quantity /i^ = (cos k — — 1 

for e I— )■ G C rather than the square root itself which is not analytic in the complex 
plan. Such equation can be rewritten 

=(^) n - -I’M ( 29 ) 

V / j=l 

where the expression of the four C-valued functions (j)j{k) can be obtained explicitly. 
Here we just need to know that G U{1), (j)j{k) ^ 4’j'{k) for j ^ j' (the gaps 

are not closing, except at 6 = n/A which has to be avoided) and (pj{k) is smooth in 
k as soon as 0 < 6 < n/2 (the bands are not degenerated to a line). Note that when 
e”’*^ = (j)j{k) one has /ifc = 0 that is p = 1: this corresponds to the gap/band edge. 

The construction of the Riemann surface consists in two steps. One is the stan¬ 
dard construction of an elliptic curve and its geometric interpretation as a compact 
surface obtained by gluing two Riemann sheets corresponding to the two branches 
± of the square root. The second step, specific to unitary systems, accounts for the 
periodic structure of the phases through the dependence in exp Tie (and its powers) 
appearing in /r^. 


4.1.2 Elliptic curves 

From now on, we consider a fixed momentum k and forget about it in the nota¬ 
tions. We start with the standard construction of elliptic curves by considering the 
Riemann surface 


"^0 — {{Z-: k) ^ 


= W(Z) = — J](Z - 


i=i 



(30) 


with the (pj satisfying the previous assumptions. This is an almost standard con¬ 
struction |3S|- First, this is a compact Riemann surface since it can be covered by 
open subset that are homeomorphic to open subset in C, with holomorphic transi¬ 
tion function on the intersections. The singularities that appear when Z —)■ 0 and 
oo are not essential because the corresponding neighborhood can also be identified 
with open subsets of C. There is a nice interpretation of this curve as a double 
covering of C, or its compactihed version |33], as we will see in the next section. 

The complex square root cannot be dehned analytically over the whole complex 
plane C, nor on its compact version, namely the Riemann sphere 5 = CUoo j35]. A 
natural branch cut on C (or S) is dehned such that W{Z) < 0. We know from Eq. 
(29) that W{Z) is real for Z G U{1) and that the sign of W changes at pi where 
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/i = 0. By construction, the (j)i are located on the equator of the Riemann sphere. 
We chose the order of the such that W{Z) < 0 for Z G [</>i,02] and Z G [03,04] 
(this notation is ambiguous in U{1) but that does not matter in the following). As a 
result W(Z) is real and positive for Z on the equator of each Riemann sphere as well, 
but outside the branch cuts. To get an analytic structure for TZq, consider two copies 
of such cut Riemann spheres and S~ and set the convention that /r = +\/W > 0 
(respectively — a/W < 0) on such a region of the Riemann sphere S'^ Hence, 

by construction of these two copies, one travels from one sphere to the other by 
crossing a branch cut, such that /r goes smoothly from positive to negative values 
|35) . The square root is now analytic over TZq (instead of S) which is a torus of genus 
1 (instead of a sphere). The Riemann surface TZq is locally homeomorphic to the 
complex plane, but not globally since the topology is different. This is the price to 
pay to have smooth functions. But this also gives a direct geometrical interpretation 
of the winding of the edge states j3T], as we will see. 


4.1.3 Punctured torus for the unitary problem 


We would like to apply the general theory discussed above for Z = exp (—ie). To 
do so we hrst need to extend the real variable e to the complex plane. The phases 
e being dehned modulo 27r, a natural extension of their domain of dehnition is the 
complex cylinder C = {z = e + iri\{e^ri) G S'! x M}, as depicted in hgure Then we 
dehne the map 


V? : 


c 

z 



(31) 


that sends this cylinder to the complex plane by preserving the circles (see Fig. |^. 
However, note that this map is not analytic since it has two essential singularities 
at 7] —>■ +00 and p —>■ —oo, mapped respectively to oo and 0 in C. This is specihc 
to unitary models where the phase is 17(l)-valued whereas such singularities do not 
appear when doing E z for real energy E of Hermitian systems |3l]. These sin¬ 
gularities will stay all along the construction and will be actually necessary. Indeed, 
the image of C by (p is the Riemann sphere S* punctured of two singular points 0 and 
oo. This surface is not compact anymore, but outside this two points the function 
is still analytic such that around the image of U{1) (the equator in this picture), 
phase terms exp(—ie) can be extended in an analytic way. 

One can now apply the general method of the previous section to construct the 
Riemann surface by replacing the two cut Riemann spheres by two punctured 
cut Riemann spheres (5*)^. The Riemann surface TZ, defined as some pull-back of 
TZo by ip, 

TZ = [{fi, z) e C X C I (/r, Z = (p{z) = e“'^) G T^o } (32) 

is obtained by gluing (iS*)"*" and (5*)~ together by the branch cuts fF(e“‘^) < 0 
which correspond to the regions of the two bands e (at hxed k) as shown in Fig. 
These bands, together with the gaps g±, constitute the equator of each punctured 
sphere as they originate from the real part of the complex variable ^ (blue circle in 
hgure 1^. All the quantities used, besides the square root, only involve polynomials 
in e“'^ and e^^ that are perfectly smooth on TZ since they are so on each copy of 
S*. The square root of /i^ is an analytic function on the Riemann surface TZ and so 
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Figure 4: Sending the complex cylinder to the complex plane, or equivalently to the 
Riemann sphere with two forbidden points (essential singularities): 0 and oo. 


are the eigenvalues of the transfer matrix. In particular, one can write A = / + /r 
that corresponds to A+ when e“^^ is in a gap and the corresponding fi is positive, 
and similarly for A_ with /i negative. In this picture A contains both square roots 
and one can go analytically from one to the other using the continuation in the 
complex numbers. Consequently the edge states of the system are now described in 
an analytic way on TZ, and the full square root is given by one single formula instead 
of two. 

Geometrically, the Riemann surface TZ is a torus with four punctured points 
(see Fig. [^. Each gap is present in two copies that generates loops along the torus. 
Indeed, the two copies of g+ generate the outer loop £+ (in green) and the two copies 
of Q- generate the inner loop £- (in blue). Importantly, due to the four singular 
points O’*", cxD’*', 0“ and oo~ that cannot be crossed, these two loops are not equivalent 
(or homotopic) as in the standard torus since they cannot be deformed one to each 
other without crossing a singular point. Hence the Riemann surface TZ keeps track 
of the two distinct gaps and g^ of the unitary problem. Finally note that such 
relative position (interior/exterior) is completely arbitrary and just depends on the 
way we draw the construction of the torus. Indeed the two conhgurations of torus 
are topologically equivalent by "twisting" the full torus. 

The construction of the Riemann surface was performed at hxed quasi-momentum 
k. When varying k, the size of the gaps changes so that one obtains a family of Rie¬ 
mann surfaces TZk- Still, the topology remains the same as these surfaces are all 
homotopic one to another (as long as the gaps do not close): one can deform all this 
family to the same torus TZ for all k. 

4.2 Winding number 

When k spans Si, an edge state may cross a gap by moving from one band to the 
other one, whereas the other edge state, located at the other boundary, crosses the 
gap in the opposite direction. On the Riemann torus TZ, these two edge states span 
one of the two loops £± and can thus be qualihed as topological: the winding of their 
pair cannot change value unless the gap closes. This dehnes a topological invariant 
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7^ = (5*)+ U (5*)- 


Figure 5: Construction of the Riemann surface TZ by gluing two punctured Riemann 
spheres (5*)"*“ and {S*)~. The branch cuts, in red, correspond to the bands. They and are 
delimited by the 4>i according to Eq. (29) and separate the two gap g- (blue) from the gap 
< 7 + (green). For convenience, we have flipped one sphere in a way such that when arriving 
at a branch cut we are going to the other surface, stay to the same hemisphere (north or 
south), as illustrated by the arrows. The gaps form two non-contractile loops and £+ 
on TZ which are non-homotopic one to each other because of the essential singularities oo^ 
and 0^. 


for each gap that counts algebraically the number of chiral edge states. As noticed 
by Hatsugai jST], this winding number is nothing but the intersection number of 
the curve spanned by the edge state dispersion relation e(fc) (after having identified 
all the Riemann surfaces TZk) and some "vertical line" of the torus, e.g. one of the 
branch cuts. 

To compute explicitly this number W, it is particularly convenient to deal with 
one continuous function e{k) for the dispersion relation of the pair of edge states 
rather than a multi-valued function as it is the case when the two edge states cross 
in the gap (see figure]^ (a)). Such a single-valued function can always be obtained by 
shifting one of the two edge state’s dispersion relation, as shown in figure [^(b). This 
is performed by adding a phase-shift to the boundary vector (e.g. bi) coefficients 
as t[ —)■ r(e*‘^ and p'l —)■ Pie'*^ so that the unitarity of Ui (and then U) is preserver^ 
Importantly, the transfer matrix is not affected by this transformation, hence both 
the bulk bands and the construction of the Riemann surface remain the same. 

It is enough to focus on the gap g- only. The single-valued dispersion relation 
e_(/c) of the pair of edge states is smooth and periodic, so that the winding number 
of exp(—ie_(fc)) always vanishes. However in the Riemann surface picture it can be 
used to wind around nontrivial loops. In this gap, e_(fc) touches the bands in 
and k 2 , that is e_(fcj[‘) = G^^t{kl) and e^i^k^) = G^^i^k^), for G^^t/b given in (14). 


That way, one gets > 0 for k E]kl, k^l and pk < 0 for k e]0, kl[U]k 2 , 7i[. Let us 


^Similarly the unitary transformation ti 
dispersion relation. 


rie“^ and pi —>■ pie“^ shifts the other edge state’s 
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Figure 6: (a) Bulk bands and edge states obtained in the thermodynamic limit with the 
transfer matrix method for 9 = vr/S. (b) The boundary state at edge re = 1 (blue) has been 
shifted ((/? = tt) so that the two edge states touch at k\ and k^- The dispersion relation of 
a pair of edge states (in each gap) is single-valued, (c) Appearance of trivial edge states 
for 9 = Sir/8 when complex phases are added to the model ([^. 


then define the function 

F_(A;) = sign(/ik)D_(k) (33) 

where D-{k) measures the relative distance between the dispersion relation of the 
edge state and a bound of the gap g^. Let us choose the lower bound as a conven- 
tiorj^ in that case, D_{k) = (e_(fc) — G_^b{k)) / {G^^t — G-fi{k)). The function F_ 
is well defined as soon as the gap does not close {6 ^ 7r/4). It is 27r-periodic in k, 
continuous in kl (since it is 0) but not continuous at k^. However the function 

re;_(fc) (34) 


is continuous, and has non trivial winding number 


= - 


wZ\k)dw_{k) = -- / dF_{k) . 


(35) 


Since FL is periodic, then, if F_ were also continuous, its winding number would be 
0. However F_ is only piecewise continuous, such that 


W\w] = - - 


fc* 


'2 1 

dF{k) - - 




dF{k) = D{kl) - D{k;) = 1 (36) 

where we have dropped the gap index because of the gap correspondence discussed 
in section 3.2.3 From an effective point of view, the computation of W[w\ amounts 
to count the number of times clfc/ifc changes sign at points k* such that F{k*) is 
discontinuous, that is D{k*) = 1. On can thus write W[w\ as an intersection number 


Wlw] = 


sign 


d/J-k 


k*s.t.D{k*)=l 


dk 


( 37 ) 


k=kf 


^Note that the other convention would reverse the sign of the winding number. 
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where the derivative dk^ik is well dehned on the Riemann surface. In practice, when 
going from a lower edge state (/ifc < 0 ) to an upper one (/ifc > 0 ) while increasing k, 
the derivative of /i*, is positive and the winding number increases by 1 , and conversely 
decreases by 1 when going from the upper to the lower Riemann sheet. 

It is clear that this analysis is still valid for several edge states and without the 
gap correspondence; the invariant W thus counting the number of topological edge 
states in a given gap. Besides, it also distinguishes topological edge states from non- 
topological "accidental" edge states that may appear for certain set of parameters, 
as illustrated in £gure[^ (c). Indeed these edge states are associated to contractile 
loops on TZ and thus do not benefit of any topological robustness. Their winding 
(36) is clearly zero. 

Finally note that the loops around the essential singularities of TZ might be 
related to previous invariant W. For example, a loop around 0“ on Figure is 
homeotopic to the one passing by —)■ 02 —t 03 —>■ 04 —>• 0 i and traveling along 
left solid red curve, dashed blue part of i~, right solid red curve and solid green 
part of £■*■, respectively. This loop corresponds to a global path in the spectrum, 
crossing both bands and both gaps but with opposite localization on the edges. In 
particular, the sum of loops around facing singularities from distinct Riemann sheet 
(e.g. 0“ and oo+) are homeotopic to U £~ and hence related to gap invariant W. 


5 Application to physical models 

As mentioned in section]^ the unitary matrix U{k) dehned in Eq. (|^ maps on 
several physical two-dimensional systems in a cylinder geometry. First, let us notice 
that from the factorized form ([^, it turns out that U actually describes an oriented 
square lattice similar to the Ho-Chalker model |36] as depicted in Fig. In that 
case, Ui and U 2 can be interpreted as scattering matrices that describe the coherent 
rehection and transmission processes at the nodes of the network. It was shown by 
Chong and collaborators that such an oriented network actually also describes the 
propagation of electromagnetic modes in arrays of optical pSl EU] or micro-wave [30] 
resonators beyond the tight-binding modeQ 

Interestingly, this model of a (static) network, also maps on other dynamical 
Floquet systems as we now show. To see it explicitly, let us factorize U{k) and 
replace the scattering parameters by their value (Eq. (|^). One gets 

U^(k) = e‘t-- , U 2 = . (38) 

Clearly, U{k) reveals a quantum protocol that consists in six steps acting on a two- 
level system. By repeating periodically this protocol, U{k) can be interpreted as 
the Floquet operator (evolution operator after one period of time) of a discrete-time 
quantum walk (see Ref. [27] for a pedagogical introduction) of a quantum system 
consisting of spin-1/2 particles located at the nodes of a lattice. The hrst hve steps 
are given by Ui{k) which is block-diagonal in the basis of the position across the 

^Up to another choice than ([^ for the scattering parameters. Despite a few technical differences 
the same method can be applied, leading to the same conclusions. 
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(a) (b) 



(d) 




Figure 7: (a-b-c) Three different physical models on a strip periodic in the x direction 
described by the same unitary matrix U{k). (a) Oriented network where the white (black) 
nodes are characterized by scattering matrices Ui (C/ 2 )- The arrows represent the incoming 
and out-going states and a phase e/2 is accumulated between two adjacent nodes. The black 
squares constitute the boundaries where a total reflection occurs, (b) Spin-1/2 particles 
on a lattice on which a protocol of six steps is applied. Five of the steps are encoded into 
Ui which is local in the transverse direction, whereas the rotation U 2 is non-local, (c) Four 
steps of a periodic hopping process. The arrows represent & 7x12 phase associated to the 
hopping, (d) Sketch of the pulses associated to each step of the cycle of period T. 


cylinder’s width. Thus, these operations are local in position as sketched in Fig. 

(b) and correspond successively to various shifts and spin-rotations. 

Finally, the non block-diagonal operator U 2 is applied, so that the corresponding 
operation (a spin-rotation by an angle 26 — n around the y axis), is applied on a 
two-level quantum state which is delocalized on sites n for spin down and n -|- 1 for 
spin up. During this last step, one spin at each edge of the strip is left unchanged. In 
contrast with the oriented network model which is static, this describes a (Floquet) 
dynamical process. 

This model can hnally be mapped onto a time-dependent tight-binding model 
on a square lattice, where the hopping amplitudes are successfully switched-on and 
off, in the spirit of previous Floquet toy models [T8l l26| |27j as illustrated in Fig. 

(c) . Unlike the two previous models, one can even write down explicitly the corre- 
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spending Bloch Hamiltonian 




Jig+(A;) • O’ 

= H^{k) 

for 

0 < t <ti 

J2^y 

= H2 

for 

ti <t <t2 

J 3 g_(A:) • a 

= Hs{k) 

for 

t2 < t < 

Jy 

= 774 

for 

ts < t < T 


(39) 


where g±{k) = (cosfc/2, ± sin/c/2, 0). By giving a step prohle to the time evolution 
of the couplings as depicted in Fig. (d), then Ui{k) and U 2 can be seen as the 
evolution operators 


dtHsik) g-i/// d.tH2 g-i/o‘^ dtHiik) 


Uo = 


■ CT 
— o ^JT-t 


dt 


(40) 


The Hamiltonian (39) has a simple interpretation sketched in Fig. (c) and (d). 
First, a coupling of amplitude Ji is switched-on between second nearest neighbors 
(A and B sites). Next, a purely imaginary intra-dimer coupling is switched-on with 
an amplitude J 2 . Then another coupling of amplitude J 3 is switched-on between 
second nearest neighbors. And hnally, a purely imaginary inter-dimer coupling is 
switched-on with an amplitude J 4 . During this last step, sites at each side of the strip 
are left uncouplecQ Considering that these four steps are repeated periodically in 
time, then the unitary matrix U{k) is nothing but the Floquet operator in a cylinder 
geometry and the phase e can thus be interpreted as the quasi-energy of the periodic 
dynamics, as for the previous model. 

These three models, described by the same unitary matrix U{k) and thus the 
same transfer matrix, share the same topological properties independently of whether 
a Hamiltonian or a periodic dynamics can be associated to them. This illustrates 
the generality of the framework we have used along the paper. 


6 Discussion 


The transfer matrix approach allows the dehnition of a gap topological invariant, as 
opposed to a bulk topological invariant dehned for the bands of the periodic system, 
e.g. the Chern number. This is particularly useful for unitary systems for which the 
Chern numbers can vanish while the system still exhibits topological edge states, as 
illustrated all along this paper. 

The topological nature of the edge states is revealed by the Riemann surface 
when taking into account the boundary conditions. Note that the construction of 
the Riemann surface only requires the transfer matrix that describes the (projected) 
bands in the thermodynamic limit. It thus only contains information about the 
bulk (one-dimensional) system. This interplay between the bulk information and 
the boundary conditions is encoded into the quantity Q{k,e) (dehned in Eq. (21)) 
whose change of sign for every e{ko) in a gap at hxed ko guaranties the existence of 
an edge state, as shown in hgure It follows that a topological transition occurs 
when the gap closes at (/cq, Cq) if and only if Q{ko, Cq) changes sign. 


^Notice that by imposing instead periodic boundary conditions, one gets the time ordered 
evolution operator for the bulk system from which a bulk gap topological invariant can be computed 
[261 [?7] and is found to be —1 (0) in the topological (trivial) phase in agreement with our approach. 
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Besides, the topological invariant W appears explicitly as an obstrnction in defin¬ 
ing a continnous and periodic argument for the phase eigenvalues e on the underlying 
Riemann surface 7Z. This invariant takes a geometrical meaning by counting the 
number of times a pair of edge states (located at two opposite edges) winds around 
one of the two loops of the punctured torus 7Z. The fact that these two loops are 
non-homotopic results from the existence of essential singularities that are absent 
in Hermitian systems for which a similar construction was first performed |31] . 

Other gap invariants have been proposed in two dimensions, especially in the 
context of Floquet systems which are described by a time ordered evolution operator 
[2S1I2Z]. In that case, the additional time parameter allows for the definition of a 
gap topological invariant for the bulk (two dimensional) system, irrespective of the 
boundary conditions. In contrast, the approach developed in this paper applies for 
any unitary system in a finite size (cylindrical) geometry, irrespective of the existence 
of a time dependent dynamics. Furthermore, other topological indexes associated 
to edge states of Hermitian systems have been already proposed |3Sl ISI- Their 
generalization to unitary systems and the relation with the present index is a natural 
direction for future investigations. Besides, it would be interesting to bridge the 
transfer matrix approach with the scattering matrix approach that also provides a 
topological characterization of edge states in both Hermitian [381 [39] and unitary 
|4U] systems. Unlike the transfer matrix which somehow probes a bulk property of 
the reduced (one-dimensional) system, the scattering matrix requires to connect the 
system to a lead and thus only probes the edges properties. Moreover, the winding 
number that is defined is assigned to one edge whereas the two edges are required 
in the Riemann torus picture. 

Finally it would be very interesting to adapt this method to other exotic topolog¬ 
ical phases, in particular when time-reversal symmetry [2ziai] or dissipation [121115] 
plays a crucial role. 
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